The thermopower and Nernst Effect in graphene in a magnetic field 
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We report measurements of the thermopower S and Nernst signal S yx in graphene in a magnetic 
field H. Both quantities show strong quantum oscillations vs. the gate voltage V g . Our measure- 
ments for Landau Levels of index n/0 are in quantitative agreement with the edge-current model 
of Girvin and Jonson (GJ). The inferred off-diagonal thermoelectric conductivity a yx comes close 
to the quantum of Amps per Kelvin. At the Dirac point (n = 0), however, the width of the peak in 
a yx is very narrow. We discuss features of the thermoelectric response at the Dirac point including 
the enhanced Nernst signal. 



In graphene, the linear dispersion of the electronic 
states near the chemical potential fi is well described 
by the Dirac Hamiltonian. As shown by Novoselov et 
al. p], 0, Q and by Zhang et al. [J, H 111 quantization of 
the electronic states into Landau Levels leads to the in- 
teger quantum Hall Effect (QHE). Because of the linear 
dispersion, the energy E n of the Landau Level (LL) of 
index n varies as E n — sgn(n)y / 2ehvpB\n\, where B is 
the magnetic induction, vp the Fermi velocity, e the elec- 
tron charge, and h is Planck's constant. The quantized 
Hall conductivity is given by 



4e^ 
~h 



(1) 



where the factor 4 reflects the degeneracy g of each LL 
(2 each from spin and valley degrees). 

Detailed investigations of the longitudinal resistance 
R xx and Hall resistance R X v have been reported by sev- 
eral groups [l|, 12, 13, |j, la, Iy, ll| • By contrast, the ther- 
moelectric tensor Sij is less investigated. Sij relates the 
observed electric field E to an applied temperature gra- 
dient -VT, viz. E =S -(-VT). On the other hand, 
the charge current density J produced by — VT is ex- 
pressed by the thermolelectric conductivity tensor a, viz. 
J = a ■{— VT). Although J is not measured directly, ay- 
may be obtained by measurements of both Sij and the 
resistivity tensor pij = Rij. (By convention, —S xx is the 
thermopower S; we refer to S yx = E y /\X7T\ as the Nernst 
signal.) 

A most unusual feature of the thermoelectric response 
of a QHE system (for n ^ 0) is that, despite the domi- 
nance of the off-diagonal (Hall- like) current response, the 
thermopower displays a large peak at each LL whereas 
the Nernst signal is small. In the geometry treated by 
Girvin and Jonson Hd (Fig. Hk), the 2D sample is of 
finite width along x, but is infinite along y (with the ap- 
plied magnetic field H||z). As we approach either edge, 
the LL energy E n rises very steeply (bold curve). At 
T = 0, edge currents I y exist at the intersections (open 
circles) of E n with the chemical potential /1. In a gradient 
— VT, the magnitude of I y is larger at the warmer edge 
than at the cooler edge because of increased occupation 
of states above y,. The difference \5I y \ is a maximum 
when [i is aligned with E n in the bulk. The correspond- 



ing value of a yx is then a universal quantum (kge/h) In 2 
with units of Amperes per Kelvin (fee is Boltzmann's 
constant). In turn, 5I y produces a quantized Hall volt- 
age Vh — (h/e 2 )SI y that drops ||x. Hence, conflating 
these 2 large off-diagonal effects, the thermopower S be- 
comes very large when [i aligns with E n (n =/= 0). By 
contrast, the transverse (Nernst) voltage is small (in the 
absence of disorder). 

In graphene, this picture needs revision when \l is 
at the Dirac point. For the n = LL, the nature 
of the edge currents is the subject of considerable de- 
bate [0, El, [H, [3. What are the profiles of S and 
Sy X ? We have measured Sij and Rij to a maximum H 
of 14 T at 20 and 50 K. Our results reveal that, at 9 T, 
the thermoelectric response in graphene already falls in 
the quantum regime at 50 K. The inferred off-diagonal 
current response a xy is a series of peaks close to the quan- 
tum value (gkBe/h) In 2 (independent of n, B and T). We 
compare our results with the caculations of Girvin and 
Jonson (GJ), and discuss features specific to the n = 
LL at the Dirac point. 

Kim and collaborators Q have pioneered a litho- 
graphic design for measuring the thermopower of carbon 
nanotubes. We have adopted their approach with minor 
modifications for graphene. Using electron-beam lithog- 
raphy, we deposited narrow gold lines which serve as a 
micro-heater to produce — VT and thermometers (Fig 
QJl). The latter are also used as current leads when R xx 
and R xy are measured. Above 10 K, the thermometers 
can resolve ST ~ ±1 mK. The typical 5T is ~10 mK 
between the 2 thermometers. However, the small spac- 
ing of the voltage leads (~2 /im) leads to an uncertainty 
of ±10% in estimating ST between them. A slowly os- 
cillating current at frequency to is applied to the heater 
and the resulting thermoelectric signals are detected at 
2ui and —90° out of phase. 

In our geometry with — VT| |x, H| |z, the charge current 
density J is (summation over repeated indices implied) 



Ji 



jEj +Oij(-djT) (i = x,y). 



(2) 



The (2D) conductivity and thermoelectric conductivity 

(e 2 /T) and oiij = 
0, we have for the 



tensors are often written 
L}j/T 2 , respectively. Setting J 
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FIG. 1: (Color online) (Panel a) The effect of — VT on the 
edge currents I y in a QHE system (n 7^ 0). The energy 
E„ of a LL (bold curve) increases very steeply at the sam- 
ple edges, with fi the chemical potential (dashed line). If 
H z > 0, I y is negative (positive) at the left (right) edge, as 
indicated by open circles. The magnitude \I y \ is larger at the 
warmer edge. Fermi-Dirac distributions f(E) are sketched 
at the sides. (Panel b) Curves of thermopower S = — S xx 
vs. gate voltage V g in Sample J10 at selected T. The curves 
are antisymmetric about the Dirac Point which occurs at the 
offset voltage Vo = 15.5 V. The peak value S m is nominally 
linear in T from 25 to 300 K (Panel c). Less complete data 
from sample K59 are also plotted. A photo of Sample J 10 
(faint polygon) is shown in Panel (d). A micro-heater as well 
as thermometers (therm) and signal leads are patterned with 
electron-beam lithography. The black scale bar is 3 /an. 
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FIG. 2: (Color online) Variation of thermopower S vs. V g 
(bold curves) and conductance G xx vs. V g (thin curves) in 
Sample J3 at H — 5, 9 and 14 T (Panels a, b and c, respec- 
tively). The offset Vo = 12.5 V. All curves were measured at 
20 K except for the curve at 50 K in Panel (c). Vertical lines 
locate the maxima of G xx . 



observed i?-fields 

E t . = -p ik akj(-djT) = Sij(-djT), (3) 

with pij = Rij the 2D resistivity tensor. The ther- 
mopower S = —E X /\VT\ equals p xx a xx + p yx a xy (S > 
for hole doping) , while the Nernst signal is given by (with 
p xx — Pyy) 

Sy X p xx & X y Pyx&xx' (4) 

The 2 terms tend to cancel mutually, except at the Dirac 
point where p yx vanishes (see below). 

Inverting Eq. [3l we may calculate the tensor ay from 
measured quantities. We have 

a xx = -(a xx E x + a xy E y )/\VT\ 

a xy = (-a xy E x + a xx E y )/\VT\. (5) 



Under held reversal (H — > -H), S is symmetric whereas 
S yx is antisymmetric. For each curve taken in held, we 
repeat the measurement with H reversed. All curves of S 
and S yx reported here have been (anti)symmetrized with 
respect to H. As for charge-inversion symmetry, we ex- 
pect the sign of S to change with the shifted gate voltage 
Vg = V g — Vo (i.e. between hole and electron hlling), but 
the sign of S yx stays unchanged (Vo is the offset volt- 
age). However, we have not imposed charge- inversion 
symmetrization constraints on the curves. Apart from 
field (anti) symmetrization, all the curves are the raw 
data. 

Figure [lb shows traces of S vs. V g at selected T. The 
thermopower S changes sign as V g crosses the charge- 
neutral point (Dirac Point) , assuming positive (negative) 
values on the hole (electron) side. The peak value S m is 
nominally T-linear from -20 K to 300 K (Fig. [Tb). 
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FIG. 3: (Color online) Variation of Nernst signal S yx vs. V g 
(bold curves) and conductance G xx vs. V g (thin curves) in J3 
at the 3 field values H = 5, 9 and 14 T (Panels a, b and c, 
respectively). All curves were measured at 20 K. Vertical lines 
locate the maxima of G xx . The sign of S yx was incorrectly 
assigned in a previous version of this paper [l5| ]. 
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In sharp contrast to the smooth variation in Fig. [TJ 
the curves of S vs. V g show pronounced oscillations when 
H is finite, reflecting Landau quantization of the Dirac 
states. Figures^, b and c display S vs. V g (bold curves) 
with H fixed at the values 5, 9 and 14 T, respectively. 
For comparison, we have also plotted (as thin curves) 
the corresponding conductance G xx = a xx . (With the 
exception of the curve at 50 K in Panel (c), all curves 
were measured at 20 K.) Whereas at large \n\, the peaks 
in S are aligned with those in G xx (vertical lines), at 
n = ±1, they disagree. In Panel a, the peaks for Vg < 
(hole doping) decrease systematically in magnitude as n 
increases from 1 to 4. 

The curves of the Nernst signal S yx are displayed in 
Fig. [3] For n ^ 0, S yx displays a dispersive profile 

centered at the vertical lines, in contrast with the peak GJ find that S displays a series of peaks, with the peak 



FIG. 4: (Color online) The thermoelectric response functions 
a xx (faint curve) and a xy (bold) vs. V g at 9 T (Panel a) and 
14 T (b) at T = 20 K calculated from S and S yx (Eq. E}. In 
(a), the width of the peak at n=0 is narrower than the others 
by a factor of 5. In both panels, vertical lines locate the peaks 
of G xx . The horizontal dashed line is (4fcse//i) ln2. 



profiles of S. Moreover, S yx is smaller in magnitude by 
a factor of 4-5. For the n = LL, however, the profiles 
change character, with S yx displaying a large positive 
peak. The reason for the enhancement of S yx at the 
Dirac point is discussed below. The positive sign of S yx 
at the n=0 LL implies that the Nernst S-field is 
parallel to H x (— VT) 

In the theory of GJ @, valid for GaAs-based 
devices, the edge current difference is SI y = 
iae/h) J2nl dk(de/dk)(df/dT)ST, with /(e) the Fermi- 
Dirac distribution. The off-diagonal term a xy is given by 

(e/hT) /~ de(e - (/,) . When (j, = E n , a xy at- 

tains a peak value, corresponding to a quantized current 
per Kelvin, given by (<7=1) 



— In 2 (~2.32nA/K). (6) 
h 



4 



value at LL n given by 

At low T, S is independent of H and T. 

In Fig. [21 the peak value of S at the n = — 1 LL in- 
creases from 25 at 5 T to 41 /iV/K at both 9 and 14 
T. Moreover, as T increases from 20 to 50 K (Panel c), 
the peak increases only weakly (41 to 48 fiV) in sharp 
contrast with the T-linear behavior at H=0 (Fig. [TJ in- 
set). We thus confirm the prediction of GJ that S (at 
the peak) saturates to a value independent of T and H 
at sufficiently low T. This saturation contrasts with the 
T-linear behavior of S in H=0 (Fig. [lb). 

In graphene, however, Berry phase effects lead to a 
the ^-integer shift in Eq. [1] [4]. In evaluating o~ xy ~ 
J2 n f(En), the ^-integer shift implies that S pea k(n) de- 
creases as ln2/(en), instead of Eq. \7\ The measured 
values Sp ea k — 41 /uV/K at n = -1 at 9 T already ex- 
ceeds slightly the predicted value 39.7 ^V/K in Eq. [7] 
Future experiments on cleaner samples may yield values 
closer to the predicted value 59.6 /iV/K for Dirac sys- 
tems. Regardless, for LL with n ^ 0, our observations 
are generally consistent with the GJ theory. In principle, 
Eq. [7] provides a way to measure ST on micron-scales 
with a resolution approaching voltage measurements. 

The most interesting question is the thermoelectric re- 
sponse of the n=0 LL. This is easier to analyze using 
the pure thermoelectric currents a xx and a xy (obtained 
using Eqs. [5]). In Fig. Uh., a xy and (a xx ) is plotted as 
bold (thin) curves for H=9 T and T = 20 K. Panel (b) 
shows the curves at 14 T. Compared with S vs. V g , the 
peaks in a xy are much narrower and clearly separated by 
intervals in which a xy is nominally zero. Likewise, the 
purely dispersive profile of a xx is also more apparent. 
Consistent with the GJ theory, the overall magnitude of 
a xy (for n ^ 0) is larger than that of a xx . 

A striking feature of a xy is that its peaks are inde- 
pendent of n. Their average value ~75 nA/K reaches to 



within 30% of the quantized value of Eq. [6] with g = 4 
(horizontal dashed line). The largest uncertainty in our 
measurement is in estimating the gradient between the 
voltage leads. For the n = LL, the shortfall may also 
reflect incipient splitting of the Landau sublevels (com- 
pare 9 and 14 T traces). The uniformity of the peaks in 
a xy accounts for the observed enhancement of the Nernst 
peak at the Dirac point. By Eq. 01 S yx is the difference 
of 2 positive terms. For the 2 terms are matched, 

and the partial cancellation leads to a dispersive profile. 
However, for n — 0, the vanishing of p yx as V' — ► 
strongly suppresses the second term —p yx a xx . The re- 
maining term p xx a xy then dictates the size and profile of 
the Nernst peak. 

We note that the sign of a xy is a direct consequence 
of the edge-currents Q. When H z > 0, I y <0 on the 
warmer edge and I y <0 on the cooler edge, as depicted in 
Fig. [TJl. As a result, a xy > for both hole- and electron- 
doping. 

Our measurements are consistent with the GJ theory 
for n t^O. For the n=0 LL, however, there is consider- 
able uncertainty about the nature of the edge states in 
graphene (or whether they exist in large H). The sim- 
ple edge-current picture for understanding the peaks in 
a xy may need significant revision for n = 0, despite the 
similarity of the peak magnitude. We also note that the 
peak at n =0 is much narrower (by a factor of 5) than 
the other peaks. This is also not understood. By going 
to cleaner samples and higher fields, we hope to exploit 
this narrow width to resolve splitting of the 4 sublevels 
at n = 0. 

We are grateful to P. A. Lee for many discussions. The 
research is supported by NSF through a MRSEC grant 
(DMR-0819860). 

Note added After we completed these experiments, we 
learned of 2 thermopower and Nernst experiments on 
graphene posted recently (Yuri M. Zuev et ai, cond-mat 
arXiv: 0812.1393 and Peng Wei et ai, cond-mat arXiv: 
0812.1411). 



[1] K. S. Novoselov et al. Science 306, 666-669 (2004). 

[2] K. S. Novoselov et al. Proc. Natl. Acad. Sci. USA 102, 

10451-10453 (2005). 
[3] K. S. Novoselov et al. Nature 438, 197-200 (2005). 
[4] Y. Zhang, J. Tan, H. L. Stormer and P. Kim, Nature 

438, 201-204 (2005). 
[5] Y. Zhang et al. Phys. Rev. Lett. 96, 136806 (2006). 
[6] Z. Jiang, Y. Zhang, H. L. Stormer and P. Kim, Phys. 

Rev. Lett. 99, 106802 (2007). 
[7] Joseph G. Checkelsky, Lu Li, and N. P. Ong, Phys. Rev. 

Lett. 100, 206801 (2008). 
[8] J. P. Small, K.M. Perez, and P. Kim, Phys. Rev. Lett. 91, 

256801 (2003). 

[9] S. M. Girvin and M. Jonson, J. Phys. C: Solid State Phys. 
15, L1147 (1982). 
[10] M. Jonson and S. M. Girvin, Phys. Rev. B 29, 1939 



(1984). 

[11] H. A. Fertig and L. Brey, Phys. Rev. Lett. 97, 116805 
(2006). 

[12] D. A. Abanin, P. A. Lee and L. S. Levitov, Phys. Rev. 

Lett. 96, 176803 (2006). 
[13] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 96, 

256602 (2006). 

[14] D. A. Abanin et al. Phys. Rev. Lett. 98, 196806 (2007). 

[15] In a previous posted version of this manuscript (cond-mat 
arXiv:0812.2866 vl), the sign of S yx was incorrectly as- 
signed. After careful checking and further measurements, 
we have determined that S yx > and a xy > at n — 
as reported here. The sign of E y is the same as that for 
the vortex-Nernst effect in superconductors. 



